Are Vortex Quasi-Crystals New Phases of Vortex Matter? 
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There seems to be a one to one correspondence between the phases of atomic and molecular 
matter (AMOM) and vortex matter (VM) in superfluids and superconductors. Crystals, liquids and 
glasses have been experimentally observed in both AMOM and VM. However, quasi-crystals also 
exist in AMOM, thus a new phase of vortex matter is proposed here: the vortex quasi-crystal. It 
is argued that vortex quasi-crystals are stabilized due to imposed quasi-periodic potentials in large 
samples or due to boundary and surface energy effects for samples of special shapes and sizes. For 
finite size samples, it is proposed that a phase transition between a vortex crystal and a vortex 
quasi-crystal occurs as a function of magnetic field and temperature as the sample size is reduced. 

PACS numbers: 74.25.Op, 74.25.-q, 74.25.Dw 



The general subject of vortex physics in superconduc- 
tors is quite interesting since there seems to be a large 
variety of possible equilibrium vortex phases in supercon- 
ductors [lj . The term "vortex matter" has been coined to 
emphasize the complexity and diversity of vortex phases 
in superconductors when compared to atomic and molec- 
ular matter. One can think of a one to one correspon- 
dence between phases in atomic and molecular matter 
(AMOM) and phases in vortex matter (VM). A liquid 
in AMOM corresponds to a vortex liquid in VM; a crys- 
talline lattice in AMOM corresponds to a vortex lattice 
in VM; an amorphous or glassy solid in AMOM cor- 
responds to an amorphous or glassy vortex system in 
VM. However, among all the possibilities discussed in the 
vortex matter literature, a very interesting one is miss- 
ing: the vortex quasi-crystal. Quasi-crystals in AMOM 
were experimentally discovered several years ago |2(| , but 
there are no corresponding experiments for vortex mat- 
ter. Thus, the present paper is dedicated to the proposal 
of vortex quasi-crystalline phases. It will be argued that 
vortex quasi-crystals maybe stabilized either by impos- 
ing quasi-periodic potentials or by boundary effects and 
surface energies in finite systems. 

The central question of this manuscript is: are vortex 
quasi-crystals new phases of vortex matter? If so, what 
suggestions can be given to experimentalists to help in 
the search for such phases? It must be emphasized that 
these questions will be addressed here on a preliminary 
basis, and further detailed work will be necessary. In or- 
der to consider the possibility of a vortex quasi-crystal 
one must seek under what conditions a quasi-crystalline 
arrangement is at all possible. A definite possibility is to 
argue that a stable vortex quasi-crystal can arise for an 
imposed quasi-periodic potential in a multilayered struc- 
ture, like for instance in the case of a superconductor 
that is grown on top of a flat quasi-crystal surface. An- 
other possibility is to create a quasi-periodic optical lat- 
tice with laser beams, trap and cool atoms, and produce 
vortex quasi-crystals in Bose or Fermi superfluids. This 
experiment is natural since the production of vortex lat- 



tices in superfluid Bose or Fermi ultra-cold atoms has 
become standard 0,0]. However, in this manuscript, we 
concentrate on the possibility of stabilizing vortex quasi- 
crystals only due to boundary effects in finite supercon- 
ducting samples, where the sample size and shape play an 
important role. This option is motivated by recent exper- 
imental studies of the vortex structure in disk, triangular, 
square, and star-shaped mesoscopic samples [HEIlIEI- 
Present technology should allow the preparation of pen- 
tagonal (pentagon-cylinders) or decagonal (decagonal- 
cylinders) samples as potential candidates to produce 5- 
fold vortex quasi-crystals. If boundaries are important 
then in the thermodynamic limit of infinite volume it 
is generally difficult to produce vortex quasi-crystalline 
order, unless a quasi-periodic potential is imposed, as 
mentioned above. 

This possibility is considered here under the follow- 
ing program. In this manuscript only the case of an 
isotropic type II superconductor in a magnetic field is 
considered. First, the bulk free energy is calculated for 
a triangular, a square and a 5-fold quasi-crystal array. 
The 5-fold quasi-crystal array is modeled by a Penrose 
tiling of the plane. It is shown that the Penrose tile ar- 
ray (vortex quasi-crystal) has a bulk free energy which 
is just a few percent higher than the triangular array. 
Second, instead of considering an infinite (bulk) system, 
a pentagon cylinder sample is discussed. The pentagon 
cylinder has a pentagonal cross-section (in the xy plane) 
with side dimension I, and with height L along the z di- 
rection. In this case, when the sample size gets smaller 
the contribution of the boundaries (surface energy) to the 
total free energy of the system becomes more important. 
The surface energy is highly sensitive to the symmetry 
and to the surface area of the boundaries. Taking into 
account the surface free energies, it is shown that the 
Penrose tile (vortex quasi-crystal) has lower free energy 
than the triangular lattice in certain regions of the mag- 
netic field versus temperature phase diagram. This is 
suggestive that a "first order phase transition" may oc- 
cur between the triangular lattice and the Penrose tiling 
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(vortex quasi-crystal) [!j . 

For simplicity, the starting point of the analysis to fol- 
low is the Ginzburg-Landau (GL) free energy density 
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for a bulk isotropic superconductor with no disorder, 
where AF S = F s — F„ is the free energy density dif- 
ference between the superconducting state (F s ) and the 
normal state (F n ), and F x = a\^{r)\ 2 + /3|*(r)| 4 /2. 
It is useful to introduce the dimensionless quantities 
/ = 9y/=0/E, p = r/A(T), A = 2<(T)A/$ , and 
TL = 2ir£(T)\(T)H/$ Q , where A is the penetration depth, 
f is the coherence length, and $o = hc/2e is the unit 
flux. Here, / = /oexp(i(/>), and A = Aa + V</()/k with 
k = A/£ is the GL parameter. Notice that when fa = 1 
the system is fully superconducting, and when /o = the 
system is normal, thus /o < 1 always. Considering the 
minimization of free energy density with respect to A and 
\l/ it is easy to arrive at equations for the dimensionless 
functions TL and / - The microscopic field is 



TL(x,y) = k 
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with H parallel to the z direction, and the equation for 
/o can be written as 



V 2 (log,?) + 2 K 2 = 0, 



(3) 



where g = /q is a positive definite function. Notice that 



TL 



K for H = H c . 2 and fo = (g = 0). The most 



general solution of Eq. ® has the form |l£ 

g(x, y) = exp [-k 2 {x 2 + y 2 )/2] exp [<y(x, y)] , (4) 

where j(x,y) satisfies Laplace's equation V 2 j(x,y) = 0. 
This means that 7(2;, y) is a harmonic function excluding 
the locations of vortices, and can be expressed as the real 
part of an analytic function of z = x + iy. This observa- 
tion has very important consequences for the microscopic 
field profile TL(x,y) of Eq. J5J), which depends strongly on 
the structure of g(x,y). 

The bulk Gibbs free energy density can be expressed 
as 11] 



G s (H,T) = G n (H,T) 
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where the parameter f3 — (g 2 ) / (g) 2 is a geometrical factor 
independent of n. The notation (• • • ) indicates average 
over volume. It is important to notice that (3 > 1 no 
matter what is the form of g(x, y) because of the Schwartz 
inequality. In addition, notice that the Gibbs free energy 
above is a minimum, whenever (3 reaches its minimum 
value. 



For the purpose of calculating the parameter (3 and the 
free energies corresponding to different vortex configura- 
tions, the analytical structure of g{x,y) in the complex 
plane is used to rewrite it as 



g(z,z) = cxp(-K 2 zz/2)|P(z)|, 



(6) 



where P(z) = NYitLii 2 ~ z %)- Here each Zi corresponds 
to a zero of g in the complex plane, and M is the number 
of zeros. The zeros Zj indicate the location of vortices. 
From now on it is assumed that there is only one vortex 
with flux $0 at each position Zj, i.e., each zero in non- 
degenerate. In this case, M corresponds to the number 
of vortices, and thus the total flux threading the sample 
is $ = A/"I>o. The normalization coefficient N just guar- 
antees that g(z,z) < 1. Depending on the locations of 
the zeros of g(z, z), it is possible to study several possibil- 
ities of periodic and quasi-periodic vortex arrangements. 
In this study only vortex crystals corresponding to trian- 
gular and square lattices and vortex quasi-crystals corre- 
sponding to the 5-fold Penrose tiling of the plane are con- 
sidered. Both the square lattices and triangular lattices 
can be generated via the tiling method, i.e., via the pe- 
riodic arrangements of identical square tiles or identical 
lozenges of internal angles (60° and 120°). The Penrose 
lattice, however, requires quasi-periodic arrangements of 
two types of tiles (lozenges), one with internal angles 36° 
and 144° and the other with internal angles 72° and 108°. 
Using the representation in Eq. ©, the values of (3 for 
the triangular, square and Penrose tiling are respectively 
fk = 1-16, /3 4 = 1.18 and /3 5 = 1.22. This immediately 
indicates that the triangular lattice has lower free energy 
than the square lattice which has lower free energy than 
the 5-fold vortex quasi-crystal (Penrose tiling) , as we ex- 
pected. However, the fact that the free energy difference 
between the triangular and 5-fold vortex quasi-crystal is 
only a few percent suggests that appropriate boundaries 
can favor 5-fold symmetry as the sample size gets smaller, 
as can appropriately imposed quasi-periodic potentials. 

In order to investigate how boundary effects can 
modify the total free energy of the system, a sam- 
ple in the shape of a pentagon cylinder of side £ and 
height L is considered. Imposing that no currents flow 
through the sample boundaries leads to the condition 
n- [— iV - 27tA/$ ] ^{x, y) = at all five side faces. The 
unit vector n points along the normal direction of each 
facet of the pentagon cylinder. The boundary condi- 
tions can be translated in terms of the harmonic function 
l{x,y) as 
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where n x and n y are the x and y components of the 
normal unit vector n at each one of the pentagon cylinder 
side faces. The solution for this boundary value problem 
can be obtained using a Schwarz-Christoffel conformal 
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map of the pentagon into a semi-infinite plane 



dz 
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where the vertices of the pentagon located at Zi are 
mapped into the points {x\, ±x%, ±£3) on the real axis of 
the w-plane. The full solution of this problem is cornpli- 
cated, and requires heavy use of numerical methods [12| • 
However, the free energy density can be estimated when 
the bulk solution in Eq. @ is treated as a variational 
solution of the boundary value problem determined by 
Eqs. ©, (J7J and (JSJ. The Gibbs free energy density dif- 
ference AG = G3 — G5 between the triangular and the 
5-fold quasi-periodic Penrose structure then becomes 



AG 
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where 0* = 0305/(05 — 0s), H c is the thermodynamic 
critical field, and e = (a^ — ct5)(2R e + + a 5 )/2R 2 ., 
with «3 = 0.93ao, 0:5 = 0.90ao, and ao = \J $o/H. The 
effective length of the sample R e — Ijyjl — t 2 , where 
t = 2cos(7r/5) is the golden mean. This expression for 
AG is valid only when H ^> §o/R 2 . The second term in 
AG takes into account the boundary mismatch energy, 
and indicates that as the size of the pentagon cylinder 
gets smaller it becomes more favorable to have a 5-fold 
quasi-crystal rather than a regular triangular lattice. No- 
tice, however, that when R e — > 00 the triangular lattice 
has lower Gibbs free energy as it must, and no transition 
to a 5-fold quasi-crystal occurs. Thus, this possible tran- 
sition may occur for finite sized samples only. From the 
condition that AG = we obtain 



H, 



1 
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where the transition 
= V2^ 
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Here, 
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to a quasi-crystal occurs. 
The phase diagram for a SU- 
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perconductor with k = 20, H C2 (0) 
R e = 10~ 6 m is shown in Fig. 
H C2 (0)[l~(T/T c ) 2 ]. 

The jump of the magnetization AM (— M3 — M 5 ) as 
a function of temperature at the critical field Hq can be 
calculated from the Gibbs free energy leading to 



AM = 

An 
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where 7* = n*y/0*e and AM S = (a 3 — a 5 )(R e + 03 + 
a 5 )/4n 2 R 2 e (l ~ 7*). A plot of AM is illustrated in Fig. 
12 for the same parameters of Fig. ^ Using these param- 
eters produces jump discontinuities AM 1=3 — 0.060G at 
T = 0, and AM ps -0.028G at T = 0.8T C . (Measure- 
ments of AM may require the preparation of an array 
of identical pentagonal cylinders to enhance the value.) 
Notice that AM < indicates that the 5-fold vortex 
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FIG. 1: H — T phase diagram for a pentagonal cylinder cut 
out of a regular cylinder of radius for R e — 10 _6 m. The solid 
and dashed lines represent H C2 and Hq, respectively, and the 
superconductor is assumed to have k = 20, H C2 (0) = 10T. 
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FIG. 2: The jump discontinuity AM = M 3 -M 5 at the critical 
field Hq for various reduced temperatures T/T c , using the 
same parameters of Fig. 



quasi-crystal is denser than the triangular vortex lattice 
at Hq, being at best a few percent denser at T = 0. 

In addition to magnetization measurements, it is also 
interesting to perform calorimetric experiments. How- 
ever, specific heat measurements are very difficult be- 
cause they require large samples. Since sample size is im- 
portant for the present discussion it is not clear that such 
experiments can be successfully performed. Neverthe- 
less, the thermodynamic relationship between the mag- 
netization and entropy jumps is revealed in the Clapey- 
ron equation AS = #3 - S 5 = -AMdH Q /dT. Since 
dH Q /dT < 0, and AM < implies that AS < 0, the en- 
tropy S3 of the triangular vortex lattice is less than the 
entropy S5 of the 5-fold vortex quasi-crystal, indicating 
that latent heat L = TAS is required to cause this phase 
transition. 

Thermodynamic quantities can provide a good under- 
standing of properties averaged over the entire sample. 
However, the use of local probes is much desired in or- 
der to reveal the change in structure from a triangular 
vortex crystal to a 5-fold vortex quasi-crystal. Thus, neu- 
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tron scattering, Bitter decoration and scanning tunneling 
microscopy (STM) experiments can help elucidate the 
structure of the vortex arrangement. 

In neutron diffraction experiments periodic or quasi- 
periodic variations ofH(x,y) will result in Bragg peaks. 
The position of these peaks determine the characteristic 
length scale of the vortex structure and its symmetry. 
The neutron scattering amplitude in the Born approxi- 
mation is 

M f 

&(q) = / ^nH{r) exp(iq • r)dr, (12) 

where fi n — 1.9leh/M n c is the neutron magnetic mo- 
ment and the M n is the neutron mass. The scattering 
amplitude 6(q) is directly proportional to the Fourier 
transform H(q) [H(q x , ?»)] of the microscopic field H(r) 
[H(x, y)} of Eq. (0. The neutron scattering cross section 
a(q x ,q y ) = 4ir 2 \b(q x , q y )\ 2 has sharp peaks at (q x ,q y ) = 
(0,0) (central peak) and at {q x ,q y ) = {±q xNm ,±q yNm ) 
(first Bragg peaks), where q X Nm = Qn cos(mir/N) and 
qyNm — Qn sin(m7r/JV), with m = 0, 1, N—l. For the 
triangular lattice N = 3 the first Bragg peak occurs at 
|Qs| = 2.31 x ir/ds, where efa is the lattice spacing. For 
the 5-fold vortex quasi-crystal (Penrose Lattice, N = 5) 
the first Bragg peak occurs at |Q 5 | = 2.46 x ir/d 5 , where 
c?5 is the side of a tile. Since the sample size is impor- 
tant for the observation of a 5-fold quasi-crystal, neutron 
scattering experiments may be difficult to perform. Thus, 
Bitter decoration or STM may be better techniques. For 
instance, STM scans at different fields and temperatures 
in the vicinity of Hq (T) should reveal the real space loca- 
tions of vortices, which can be Fourier transformed (FT) 
to obtain a 6-fold pattern for the triangular vortex lattice 
and a 10-fold pattern for the 5-fold vortex quasi-crystal 
(Penrose lattice), see Fig- El 




FIG. 3: First peaks of the square of the FT pattern for the 5- 
fold vortex quasi-crystal (Penrose lattice). Notice the 10-fold 
symmetry revealed. 

Now that the phase diagram, thermodynamics, and 
the local signatures of a vortex quasi-crystal have been 
discussed, it is important to say a word on the stability 
of such structures. A stability analysis in the free en- 
ergy can be performed by moving the vortices away from 



their equilibrium positions z% to Zi + 5zi. The eigenval- 
ues associated with these displacements indicates that 
for H > Hq(T) the vortex quasi-crystal lattice is sta- 
ble. However, there is no full compatibility of the 5-fold 
vortex lattice with the pentagon cylinder geometry, and 
thus the appearance of disclinations and dislocations is 
possible. As a result there is an additional possibility of 
the coexistence of a solid (crystal or quasi-crystal) vor- 
tex phase in the center of the sample and a liquid vortex 
phase closer to the edges. 

In conclusion, we have shown that vortex quasi-crystals 
may be experimentally observed in isotropic type II su- 
perconductors, provided that the sample size and shape 
are properly chosen. This result opens the possibility of 
a new phase of vortex matter: the vortex quasi- crystal. 
This yet non-observed phase adds to the many phases 
of vortex matter already seen experimentally, the vortex 
crystal, the vortex liquid [pHHIi) and the vortex glass (lif . 
By taking into account boundary effects, sample shape 
and size or properly imposed quasi-periodic potentials, 
the present work suggests that a vortex quasi-crystal 
phase may exist and compete with vortex liquid, glass or 
crystal phases, as magnetic field, temperature and disor- 
der are varied. 
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